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td : ABSTRACT 

^ ■ A previous paper proposed a new kind of random walk on a spherically-symmetric 

►-^ 

, , , lattice in arbitrary noninteger dimension D. Such a lattice avoids the problems associated 

rS 
d ' with a hypercubic lattice in noninteger dimension. This paper examines the nature of 

spherically-symmetric random walks in detail. We perform a large-time asymptotic analy- 
sis of these random walks and use the results to determine the Hausdorff dimension of the 
process. We obtain exact results in terms of Hurwitz functions (incomplete zeta functions) 
for the probability of a walker going from one region of the spherical lattice to another. 
Finally, we show that the probability that the paths of K independent random walkers 
will intersect vanishes in the continuum limit il D > ^_i - 



PACS numbers: 05.40.+J, 03.29+i, 02.50+5 



I. INTRODUCTION 

In a recent paper ^ we devised a D-dimensional spherical lattice on which a random 
walk is well defined even when D is noninteger. The advantage of such a lattice is that 
for all real values of D the probabilities have sensible and acceptable values; that is, they 
lie between and 1. A hypercubic lattice gives acceptable probabilities only for integer 
D. Having constructed such a model, our objective in the current paper is to examine the 
nature of spherical random walks in much greater detail. 

The simplest spherically-symmetric geometry is a set of regions bounded by concentric 
nested spheres. Such a geometry does not single out the integer dimensions as being special 
in any way. Thus, a random walk can be defined in any real dimension D. We consider 
an infinite set of concentric nested spheres of radius Rn {n = 1, 2, 3, . . .). Region n is 
the volume between the n — 1st and the nth spheres (see Fig. 5 of Ref. 1). At each time 
step a random walker in region n must move inward or outward to an adjacent region; the 
probability of walking outward is denoted -Pout(^) and the probability of walking inward is 
Pin(n). In our model the probabilities -Pout(^) and Pin{n) are determined by the relative 
surface area between adjacent regions. Thus, 

P,^,{n) = ^DJ^rr) 

SoiRn) + SniRn-l) ,. .N 

PUn) = SpiRn-i) 

SoiRn) + SD{Rn-l) 

where Sd{R) is the surface area of a L)-dimensional sphere of radius R: 

Since the surface area Sd{R) is proportional to R^~^, the inward- walk and outward- walk 
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probabilities are 



Poudn) ^^ 



D-l 



D-1 , r>D-l ' 



ttn -\- it _ 1 

.Z.-1 (1-3) 

p. (n) = l^!i^^ 

These formulas are valid for all D and all n > 2. For n — 1, Pout(l) = 1 and -Pin(l) = 
because there is no region inward from n = 1 and therefore the walker must go outward. 
Note that Pout and Pin are nonnegative numbers less than or equal to 1 with Pin(n.) + 
Pout(?^) = 1; thus, they are acceptable probabilities for all real D. 

A spherical random walk is described by the function Cn,t, which represents the prob- 
ability that a random walker arrives at region n at time t. The probability Cn,t satisfies 
the recursion relation 

C^,t = Pin{n + l)Cn+i,t-i + Pont{n - l)Cn-i,t-i {n > 2) , 

(1.4) 

C'l,t — Pin(2)C2,t-l , 

where Pin and Pout are given by (1.3). This equation states that for a random walker to 
be in region n at time t, she must have come from region n + 1 or region n — 1 at time 
t-1. 

In Ref. 1 we considered the consequences of the particular initial condition that at 
time t = the random walker is located in the central sphere {region 1). We represent 
this initial condition by 

Cn,0 = Sn,l ■ (1-5) 

Two immediate consequences of (1.4) and (1.5) are that for all values of t, 

oo 
I]C'n,t = l (1.6) 
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and that Cn^t > for all n > 1 and t > 0. Hence, for all n and t, C^^t lies between and 
1 and is thus an acceptable probability for all real D. 

In general, for arbitrary D the solution Cn,t to the partial difference equation in 
(1.4) and (1.5) is rather complicated. However, in Ref. 1 we used the theory of continued 
fractions to solve (1.4) and (1.5) for Ci,2t for any value of D and for any sequence of radii 
Rn- We simply iterate (1.4) and (1.5) to obtain 

Ci,o = 1 , 

C'i,2 = Qi , 

(1.7) 

Ci,4 = gi(Qi + g2) , 



Ci,6 = Qi[(Qi + Q2)' + Q2Qs 



and so on, where 



Q^ = P,^{n + l)Po^t{n) . (1.8) 

The sequence of formulas in (1.7) is a known pattern:^ Ci,2t are the coefficients in the 
Taylor series for the function whose continued-fraction representation has the coefficients 

oo 

J2 Ci,2tx' = 1/(1 - Qix/{1 - Q2x/{1 - Qsx/{1 - Q^x/ . . .)))) . (1.9) 

Thus, at a; = 1 we have the formal result 

oo 

J2 Ci,2t = 1/(1 - Qi/(i - Q2/(i - Qs/{i -Qa/.. .)))) • (1-10) 

If we truncate the continued fraction in (1.10) after Qn and substitute 

>D-1 



O - ""^ 

ri^ -\- -rt2 

^2D-2 

{R^_l + Rn ){Rn + Rn+l) 
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(1.11) 



"■ / 7-> \ D-1 



an amazing algebraic simplification occurs for each value of n: 

1/(1 -Q^/{l-Q2/{l-Qs/{l-Q4/... {I -Qn)..m = J2{^y ' ■ (1-12) 
We may then take the limit n -^ oo if the series in (1.12) converges. We obtain 

OO °° / E> \ D — 1 

For the special case Rn = n, corresponding to equally-spaced concentric spheres"^, and 
D>2, (1.13) becomes 

OO 

j2ci,2t = aD-i) , (1.14) 

where ( is the Riemann zeta function. Note that the sum in (1.14) is not a probability 
because the probabilities Ci,2t are not disjoint. Nevertheless, the sum does have a physical 
interpretation. It is the mean time spent by the random walker in region 1.^ 

Another interesting probability function is P2t {t = 1, 2, 3, . . .), the probability that 
the random walker returns to the central region (region 1) at time 2t for the first time. 
The generating-function formula^ 



OO ^ 

J2P2tx' = 1- OO ^ 



2tX* 



states the relationship between Ci^2t and P2t- The generating-function formula in (1.15) 
is universal; it holds for a random walk in any geometry. 

If we sum over all disjoint (mutually exclusive) probabilities P2t we obtain Hd, the 
probability that a random walker will eventually return to her starting point in region 1: 



UD=J2P2t . (1.16) 



t=l 
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To compute Hd we set x = 1 in (1.15) and use (1.10) and (1.13): 

Ud = Qi/(i - Q2/(i - Qs/{i - Q4/(i -Qb/.. •)))) 

_ _ 1 (1.17) 

Y^oo Ri \ 

2^N=1 \Rr,J 

Thus, for the special case in which the random walk occurs on a geometry of concentric 

spheres whose radii are consecutive integers, Rn = n, we have 

Ud = 1 {D<2) , 

1 (1-18) 

CiD-l) "- ^ 

In this paper we extend the analysis of the partial difference equation (1.4). In Sec. II 

we perform a detailed asymptotic analysis of (1.4) for large t. We use the results of this 

analysis in Sec. Ill to compute the spatial moments of Cn,t- The spatial moments represent 

the random walker's mean distance from the origin at time t. From the spatial moments 

we can show that the Hausdorff dimension of a spherical random walk is 2 for all -D > 0. 

In Sec. IV we show how to calculate the temporal moments. This calculation allows us to 

determine the average time it takes for a random walker to return to the central region. 

In Sec. V we generalize Ud to a quantity H^'", the probability that a random walker who 

begins in region n will eventually reach region m {n > rn). We show that H^'" can be 

expressed in closed form in terms of a Hurwitz function (an incomplete zeta function). 

Finally, in Sec. VI we examine the problem of intersecting random walks. We show that 

in the continuum limit the probability that K independent random walks will intersect 

vanishes when the dimension D exceeds 2K/{K — 1). This reproduces the value of the 

critical dimension for which a scalar cf)'^^ quantum field theory becomes free. 



II. LARGE-TIME ASYMPTOTIC ANALYSIS OF SPHERICAL RANDOM 
WALKS 

The partial difference equation (1.4) augmented by the initial condition (1.5) deter- 
mines the function Cn,t that describes a random walk beginning in the central region of a 
spherically-symmetric lattice. In this section we perform a large-time asymptotic analysis 
of this function Cn,t for the fiat-space case Rn = n. Substituting R^ = n into (1.4) gives 
the partial difference equation 

r c I (n-i)^-^ 

2D-1 

1 

1 -f- 2"^-i 

(2.1) 



C*!,* — -, , ^r)_i ^2,t-l 



For arbitrary D these equations do not have a closed-form solution. However, we 
showed in Ref. 1 that for D = 1 and D = 2 there are simple, exact expressions for Cn,t- 
For D = 1, 






Ci, 



2t 



(n + 2j-l)! 
(n + j-l)!j!2^+2j-2 

(2t)! 
t!t!22* 



(2.2) 



For D = 2, 



_ (2n-l)(n + 2j-l)! 
^-'-+2^-1- ^•!(2n + 2i-l)!!2^- ' ^^'^^ 

from which we have 

Ci 2t = ^^ • (2.4) 

' l + 2t ^ ' 

We begin by examining the large-time behavior of C^^t &!■ the dimensions D — 1 and 

D = 2 using the above exact solutions. We can expand the expressions in (2.2) and (2.3) 
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for large t using the following asymptotic relation satisfied by the Gamma function, 



r(i + a) 



t' 



1 + 



O? — OL 



2t 



+ o(t- 



(t ^ oo) 



(2.5) 



which may be derived from the Stirling expansion.^ We obtain for D = 1 

23/2 



Cn,t ~ 

and for D = 2 



2"^-^" + ^0(i-) 



(t ^oo) 



y/nt 






4t 




1 


1- 


1 

' 8t 


+ o(t- 


-') 



(t ^ oo, n>2) 



(2.6) 



a 



n,t 



2n-l 
t 



1- 



n — n + 1 
2i 



+ o(t- 



it 



oo 



(2.7) 



Note that these asymptotic behaviors appear to have the form of Frobenius^ series 
about t = oo. That is, they have the form of Taylor series in inverse powers of t multiplied 
by a fractional power of t. This suggests that we might seek an expression for Cn,t of the 
form 

oo 

C^,t = rC{n)Y,Ck{n)t-'' , (2.8) 

fc=0 

where we choose co{n) = 1. Actually, detailed numerical calculations for Ci^2t out to 
2t = 200, 000 show that the formula in (2.8) is not quite correct for arbitrary D; there 
appear to be logarithmic corrections which become negligible when n is large compared 
with 1. In the following analysis the parameter n will be sufficiently large that we may 
ignore these logarithmic terms. 

Frobenius series are particularly appropriate for the local analysis of differential equa- 
tions because the derivative of t°'~^ with respect to t gives another term of the same form 



with the index k shifted by one. However, they are not suitable for difference equations 



because the discrete difference in the t variable of the expression t° ^, 

(t)"-^-(t-l)"-^ , 

is not simply expressible as a combination of terms of the form t"~'^. As is shown in Ref. 5, 
an appropriate form for a Frobenius series valid in the vicinity of t = oo and useful for 
difference equations is 

Cr.,t = C{n)f^c,{n)^^^,^^ , (2.9) 

k=o ^ ' 

where co(n) = 1. Note that we have replaced t"~^ in (2.7) by r(t + a — k)/r{t). This 
ansatz has two advantages. First, for large t, we have 

r{t + a-k)/r{t) -r-'^ , (2.10) 

by virtue of Eq. (2.5). Second, the discrete difference of r(t + a — k)/T{t) is simply 
expressible in terms of the same structure with k shifted by one: 

V{t + a-k) V{t-l + a-k) . ^.V{t + a-k-l) ,^^^, 
^ = [a — k) . 2.ii 

r(t) r(t-i) ^ ^ r(t) ^ ^ 

Indeed, (2.11) is the discrete analog of the derivative relation 

^^-^ = (a - k)r-^-^ . 

dt ^ ' 

Our program now is to substitute (2.9) into the partial difference equation (2.1). To 
do so we make use of the identity 

a,_i=C(n)f;c.(n)^^^^f^ + C(n)f;(fc-a-l)c._i(n)^^^^f^ . (2.12) 
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Then we compare coefficients of r(t + a — k)/T{t). This comparison gives a recursion 
relation to be solved for C{n) and Ck{n) (n > 1). There are three cases to consider, D > 1, 
< L* < 1, and L* < 0. 
Case 1: D > 1 

This is the simplest case because the three equations in (2.1) reduce to a single equation 
valid for all n > 1; there are no special equations for n = 1 and n = 2. The /c = equation 

determines the coefficient C{'n): 

^D-i (n_i)D-i 

C(") = („ + i).-.+„.-i C(" + 1) + („_i)Li4_2)°-. C(" - 1) • (2.13) 

The solution to this equation is 

C{n) = A[n^-^ + {n - 1)^-^] , (2.14) 

where A is a constant that cannot be determined by (2.13) because the equation is homo- 
geneous. 

The k > 1 equations yield a recursion formula that gives the coefficient Ck (n) in terms 

of Cfc_i(n): 

[n^-i + {n- l)^-^]cfc(n) - n^'-^Ckin + 1) - {n - l)''-^Ck{n - 1) 

= {k-a-l)[n^-^Ck-i{n + l) + {n-l)^-^Ck-i{n-l) . (2.15) 

From the first of these equations, /c = 1, we obtain a formula for ci{n): 

n—l ^ m 

c,{n) = ci(l) + a 5^ ^^ 5^[j^-i + {3 - 1)^-^] . (2.16) 

Note that the solution for ci(n) contains the arbitrary summation constant ci(l) and the 
(as yet undetermined) index a. 
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The remaining equations express Ck{n) as a sum over Ck-i'- 

n—l ^ m 

Cfc(n) =Cfc(l) + (a + l-A;) ^ ^^ ^[j^-icfc_i(j + l) + (j-l)^-icfc_i(j-l)] . (2.17) 

lib 

For each value of k we have a new arbitrary constant; we expect an infinite number of 
arbitrary constants because we are solving a partial difference equation. In principle, 
these constants can be determined by a boundary condition. However, no such boundary 
condition is available to us because we are performing a local asymptotic analysis at t = oo. 
Ultimately, we will be interested in the large-n behavior of Cfc(n) and thus these arbitrary 
constants will prove to be irrelevant. 

We will now determine the value of the index a. We do so by examining the behavior 
of Ck{n) for large n. We already have 

co(n) = 1 . (2.18a) 

Next, we evaluate ci(n) in (2.16) for large n. The upper regions of the sums dominate in 
this asymptotic approximation and the result is independent of the constant ci(l): 

ci(n) ~ -^n^ {n» 1) . (2.186) 

Using the result in (2.18b) we find the large-n behavior of C2{n) from (2.16); again, the 
upper regions of the sums dominate and the result is independent of the constant C2(l): 

C2in) ~ °^(°^ ~ ^) !L (n » 1) . (2.18c) 

Iterating this process we find that 

, , a(a — 1)(q; — 2) n^ , . . ,x 
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In general we obtain the result 



, , a(a — l)(a — 2) ■ ■ ■ (a — k + 1) -n?^ , , , , 

Equations (2.18) can be summarized as 

If we set 

a = -D/2 (2.20) 

then (2.19) greatly simplifies: 

Ck{n) - ^ ^^^, (n»l) . (2.21) 

Note that in the above asymptotic analysis the variable n cannot exceed 0{^/t). This is 
because we are performing an asymptotic expansion valid for large t. Subsequent terms in 
the expansion (2.9) must be of decreasing size for our analysis to be correct. 

We will now argue that (2.20) is in fact the correct value of the index a. We substitute 
(2.10), (2.14), and (2.21) into the discrete Frobenius series (2.9) and obtain for large t, 

oo ( 2\k 

Cn,-'2An^-'Y.^-'-'"'h^ [t ^ oo, 1 « n < 0(v^)] , 

fc=0 

which sums to 

Cn,t ~ 2A?i^-it-^/2e-"'/(2t) (^ ^ oo, 1 « n) , (2.22) 

where the upper asymptotic limit on n is no longer necessary. Observe that as a function 
of the variable n, C^^t rises until it reaches a maximum at 

n = ^{D-l)t (2.23) 

12 



and then decays to zero exponentially. To illustrate this maximum we have computed Cn,t 
numerically for D = 3 from (1.4-5) and graphed the results in Fig. 1 as a three-dimensional 
plot. Observe that the crest follows a parabola. Note also that the maximum is very broad. 
We now impose the normalization condition in (1.6), which states that at each time t 
the probability that the random walker is somewhere on the lattice is unity. Substituting 
(2.22) into (1.6) gives 

oo 

l = 5^2A(2n)^-it-^/2e-2-'/* , (2.24) 

where we have assumed that t is odd and thus we must sum over even values of n only. 
We could just as well take t to be even and to sum over odd n only.^ Since the dominant 



contribution to this sum comes from large values of n near a/(-D — l)t we need not be 
concerned about Cn,t for small n and we may immediately replace the sum in (2.24) by 
an integral over n: 

POO 

1= / rfn2A(2n)^-it-^/2e-2-Vt . 

Jn=0 

We evaluate this integral and obtain a value for the constant A: 

2 

^ - 2^/2r(D/2) • ^^-^^^ 

It is crucial that this constant be independent of the value of t; we find that in evaluating 
the integral over n the variable t scales out and all t dependence does indeed disappear. The 
requirement that the result be independent of t is precisely the condition that determines 
the value of the index a in (2.20). Using the value of A in (2.25) and taking just the first 
term in the series (2.9) we obtain the leading asymptotic behavior of Cn^t ^or large t: 

Cn,t-^jjj^^^^^[n''-' + {n-l)^-']t-^/' (t^oo) . (2.26) 
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The result in (2.26) agrees exactly with the leading asymptotic behavior in (2.7) for 
Cn,t 8it D = 2. It also reproduces exactly the leading asymptotic behavior in (2.6) for Cn,t 
at D = 1. (To obtain the correct result at n = 1 we must be careful to set n — 1 first 
and then to allow D to approach 1 from above.) We have verified the behavior in (2.26) 
numerically by computing Cn,t for various values of D out to t = 500 and we find that our 
asymptotic approximation agrees with the exact result to many decimal places. 
Case 2: < L* < 1 

The analysis of Case 2 is more difficult than that of Case 1. If we formally assume a 
Frobenius series (2.9) and substitute this series into the partial difference equation (2.1), 
then as in Case 1, we can determine the large-n behavior of the series coefficients. We 
substitute the Frobenius series into the normalization condition (1.6) and, as with Case 1, 
we convert the sum in (1.6) to an integral that is independent of t provided that the index 
satisfies a = —D/2, the same result as that in (2.20). 

When D < 1 we encounter a problem not present in the analysis of Case 1. Recall 
from equation (2.23) that when D > 1 the function CrL,t increases with n until it reaches 



a crest at n = ^/{D^^T)t and then it decays. This crest follows a parabolic curve in the 
(n, t)-plane. The crest structure justifies replacing a sum in (2.24) by an integral, which 
we then evaluate to determine the constant A exactly. In Case 2, this crest still exists 
but it lies along a straight line near n = 1. Thus, the large-n asymptotic behavior of the 
Frobenius series coefficients is not sufficient to determine the value of the constant A. To 
illustrate the change in the crest structure we have calculated Cn^t numerically for D — 1/2 
and plotted the results in Fig. 2. Observe that unlike the situation shown in Fig. 1, the 
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crest is parallel to the t axis. 

Although we cannot derive the coefficient of the leading asymptotic behavior of Cn,t 
for large t we will make the assumption that we can analytically continue the result in 
(2.25) from values of D greater than 1 (Case 1) to values of D less than 1 (Case 2). Of 
course, the analytic continuation of an asymptotic approximation is not always the same 
as the asymptotic approximation to the analytic continuation of a function. Nevertheless, 
it is worth examining numerically the possibility that the leading asymptotic behavior of 
Cn,t for large t is given by (2.26) (supplemented by the special condition that the behavior 
of Ci^t is obtained by first setting n = 1 and then continuing D below 1). To test this 
hypothesis we have computed Cn,t for all values of t up to 3000 for D = 1/2. In Fig. 3 
we compare the exact value of Ci^t with the leading asymptotic approximation to this 
function taken from (2.26): 



93/4 
C^^^^^t-^/^ . (2.27) 

'* r 1/4 ^ ^ 



Clearly, (2.27) is an extremely good approximation for large t. 
Special Case: D = 

For this special case there is no simple formula for Cn,t'i indeed, the numbers Ci 2t are 
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quite complicated looking: 



C'l.o = 


1 , 


Cl,2 = 


2 

3 ' 


C'1,4 = 


26 

45 ' 


Cl,6 = 


502 
945 ' 


Ci,8 = 


7102 
14175 




44834 


ClAO = 


r^f^r r r 



93555 
Nevertheless, we have found an exact formula for P2t that is extremely simple:^ 

P2t = r- . (2.28) 

(2t + l)(2t-l) ^ ' 

Although it is probably valid to continue (2.26) to values of D less than 1, this con- 
tinuation clearly cannot be valid at D — because r(D/2) becomes divergent. Note that 
as D approaches the algebraic behavior in the variable t disappears; also, the coefficient 
in (2.26) vanishes because r(0) = 00. This suggests to us that at D = 0, Cn,t continues to 
vanish in the limit as t ^ 00 but that it vanishes less rapidly than algebraically. In fact, 
Cn,t exhibits a logarithmic decay as t ^ 00. 

It is easy to see why there is a logarithmic decay in Cn,t for large t. If we assume that 

C„,t ~ J^{t)C{n) {t ^ 00) 

and we substitute this behavior into (2.1) we obtain a formula for C{n) that is the exact 
analog of (2.14): 

C{n) =(- + -^] C(l) (n > 1) . (2.29) 

\n n — 1 J 
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To find C(l) we use the normalization condition in (1.6). Note that the sum over n in 
this case would be infinite except for the fact that Cn,t must vanish for n > t + 1. Thus, 
summing over n with an upper limit of t gives the crude estimate that J-'{t) ~ (lnt)~^ as 
t ^ oo. This analysis is not delicate enough to give the overall multiplicative constant in 
the asymptotic behavior. 

A more precise result can be obtained by solving (1.15) for Yl'^QX^Ci^2t and substi- 
tuting the exact formula for P2t in (2.28). We obtain 

Next we multiply (2.30) by x~^~^ and integrate around a contour enclosing the origin in 
the x-plane. We can evaluate this contour integral using the method of steepest descents.^ 
A saddle point occurs at x = 1. We find that the large-t asymptotic behavior is 



Ci,2t~^^ (t^oo) . (2.31) 



We have verified this result numerically for values of t up to 6000. 
Case 3: D <0 

When D < we find an entirely new behavior for Cn,t for large t; namely, that Cn,t 
approaches a nonvanishing steady-state function of n: 



lim Cn,t = Cn,oo . (2.32) 

t — >C!0 



To determine the function Cn,oo we merely set t = oo in (2.1) and solve the steady-state 

17 



equation: 

(n + 1)^ ^ + n-^ ^ (n — 1)^ ^ + (n — 2)-^ J^ 

2-D-i 

(-^2,cx) = oD-l I 2-D-l 3,00 +^l,cx) 7 

Ci,oo = ^_^ 21^-1 ^2,00 • 

(2.33) 

The solution to (2.33) is uniquely determined up to a single multiplicative constant Ci^oo'- 

Cn,oo = [n^-' + (?i-l)''-']Ci,oo (n>2) . (2.34) 

To find the value of Ci^oo we normalize the solution in (2.34) using (1.6) in the limit 
t -^ 00. Recall that the solution Cn^t is nonzero in a checkerboard pattern.^ Thus, we can 
take t to be odd, in which case (1.6) becomes 

00 

/ ^ C2n,t = 1 , 
n=l 

or we may take t to be even, in which case 

CXD 

/ ^ C2n+l,t = 1 • 



n=0 



In either of these two cases, after we substitute (2.34) and let t -^ 00, the normalization 
condition becomes 



(1^-i + 2^-^ + 3^-^ + 4^-^ + ■ ■ OCi.oo = 1 



whence 



^1,00 = ^^j^ {D<Q) . (2.35) 
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Thus, for large t we have^ 

Cn,t--7^^[n''-' + {n-l)''-'] (t^oo, n>2), 

^ ^ ' (2.36) 

To illustrate the result in (2.36) we present a three-dimensional plot of the function 
Cn,t in Fig. 4 for D = —1. Observe that as t ^ oo, Cn,t approaches a time-independent 
distribution. 

An alternative way to obtain the value of Ci oo rnakes use of a crucial symmetry 
property of the partial difference equation (1.4) referred to in Ref. 1 as the reciprocity 
property of (1.4). It was noted in Ref. 1 that if the value of D is reflected about the 
point D = 1 then Pout{^) and Pin{n) in (1.3) interchange roles when n > 2. From this 
observation it was shown from the properties of continued fractions^ that 



B, 



2t 



D 



(Ci,2t-2 - Ci,2t) (t>l) , (2.37) 



2-D 



where P2t is defined in (1.15). An immediate consequence of this property can be obtained 
by summing (2.37) from t = 1 to oo: 

oo 
t=l 

= (C'l.o ~ Ci 
-1-Ci 



D 



(2.38) 

2-D 



2-D 

Hence, n^ = 1 if and only if linii^oo Ci^2t = in (2 — D) dimensions. Thus, the fact 
that in Cases 1 and 2 we have Cn,t ^ as t — ;► cxd when D > implies that Ud = 1 for 
D < 2. However, we have computed Hd explicitly [see (1.17)] for D > 2. Thus, using 
(2.38) we reproduce the value of now Ci^oo given in (2.35). 
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To obtain higher-order corrections to the leading asymptotic behavior in (2.36) we 
assume that the difference between Cn^t and its asymptotic behavior in (2.36) has the 
form of the discrete Frobenius series on the right side of (2.9). (This procedure wiU be 
successful because, as we have verified numerically, the difference between (2.36) and the 
exact result has a crest as in Case 1.) We proceed exactly as in Case 1. Again we find 
that C{n) is determined up to an unknown constant A: 

C{n)=A[n^-^ + {n-l)^-^] (n > 2) , 

(2.39) 
C{l)=A . 

Note that this equation is slightly different in form from that in (2.14) because D < 0. 

Next we find a recursive formula for the coefficients Ck{n): 



Ck{n) = Ckil) - (a + 1 - /c) ^ Ck-i{r 



n-l 

(2.40) 

n— 1 m ^ ' 

+ (a + 1 - fc) 5^ m^-^ 5^ j^-i[c,_i(j) + cu-i{3 + 1)] • 

7Tl=l i = l 

This formula allows us to determine Ck{n) for large n. From our assumption that 

co{n) = l , (2.41a) 



we have 



and so on. 



ci{n)^,^ai-D)n'-^ {n»l) , (2.416) 

^^("^^ (2-1)K4-D) ^(^-^^^ ^^»^^ ^ (2.41c) 

C3(n)^, 2a(a-l)(a-2) (^_^^n^ ^^^^^^ ^2.41rf^ 
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We will now show that the choice of index 

a = D/2 , (2.42) 

which we have found from numerical computations, is consistent with the normalization 
condition (1.6). [When D < the normalization condition does not uniquely specify the 
value of a but we can at least show that the choice of a in (2.42) does not contradict the 
normalization condition.] If we substitute this value of a into (2.41), we find that 

^'("' - (2k D)(-2f-Hk-n <*-'■" >" '' ■ <'■''' 

[Note that in this equation the variable n should not exceed O(v^). See (2.21).] Thus, for 
large t the first correction to C„,oo5 apart from the multiplicative function C{n), is 

oo f-k+D/2^2k-D 

t-l^ + DCd - D) E p,_'^)(_,;-.(,_i), . (2.44) 

where we have substituted the asymptotic result in (2.43) into the Frobenius series in (2.9) 
and taken t large. The sum in (2.44) can be represented as an integral 



Jo 



(2.45) 



We now multiply the result in (2.45) by C{n) in (2.39) and sum over all even n (or 
all odd n). It is justifiable to replace the sum over n by an integral for each fixed x in the 
integral in (2.45) because there is a crest for large n. The result of doing the n integration 



IS 



At^/2^(1 -D)(l + D f dx x-^-A = 



(2.46) 



This shows that there are no small time-dependent corrections to the normalization con- 
dition in (1.6). 
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Summary of the large-time Asymptotic Behavior of C„ t 

In this section we have investigated the asymptotic behavior of Cn,t as t ^ oo. We 
can summarize this behavior as follows: 

When D > 0, we find that a = —D/2. Furthermore, for D > 1, we have determined the 
coefficient of t": 

When D < 0, we find that a — and thus Cn,t approaches a time-independent function of 
n as t ^ oo; this function is given in (2.36). Moreover, we have argued that the difference 
between Cn,t and its limiting value at t = oo is a function that falls off with increasing t 
like t^/^. We have verified this last result using detailed computer calculations. 
The special case D = has a logarithmic behavior for large t: 

/I 1 \ 2 

Cu,t ~ - + — (t ^ oo) . 

\n n — 1/lnt 

From the large-t asymptotic behavior of Ci^2t we can determine the large-t asymptotic 
behavior of P2t defined in (1.15) using the reciprocity relation in (2.37). This asymptotic 
behavior takes the general form 

P2t--Vt^ (t^oo) . (2.47) 

For all values of D 7^ 2 we derive from (2.37) that 

^ = _1_L^_^ . (2.48) 
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This result is definitely not true when D = 2 because 112 exists; it is clear that there 
must be logarithmic corrections in the asymptotic behavior of P2t for this special value of 
D. Indeed from (2.31) we see that 



2i 



(t ^ oo) . 



t[ln(2t)]^ 
When D > 2 the coefficient V is not known. However, when D < 2 we have 



- 7 



(2-L>)2^-2 

^= r(i-D/2) ■ <'■''' 

The coefficient in (2.49) can be verified for the case D = 0. From (2.28) we see that 
for large t, P2t ~ ^~^/2, which agrees with (2.47-49). Equations (2.47-49) are also valid at 
D = 1. It was shown in Ref. 1 that at D = 1 

P2t = jrTnPL^ ■ (2.50) 



For large t (2.50) becomes 



t!t!(2t-l)22* 



P2t - , As/9 ' (2-51) 



2v/^t3/2 



which agrees with (2.47-49). 
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III. SPATIAL MOMENTS OF RANDOM WALKS 

In the previous section we calculated the large-t asymptotic behavior of the solution 
to the partial difference equation (2.1) that describes a random walk on a set of concentric 
spherical shells of radii Rn = n. This random walk is subject to the initial condition (1.5) 
that states that the random walker begins in the central region. In this section we use the 
results of this asymptotic analysis to determine the mean distance of the random walker 
from the central region as a function of time t for large t. 

The kth spatial moment of a random walk is defined as a weighted average over the 
probabilities Cn,t'- 

CO 

{R\ = Y^n^Cr., . (3.1) 

n=l 

To evaluate the sum in (3.1) we substitute the asymptotic behavior of Cn,t^ assume 
that t is odd and sum over even n (or equivalently we may assume that t is even and 
sum over odd n). The leading large-t asymptotic behavior of Cn,t is completely known for 
D>1 [see (2.26)]. We thus have 

A + -D/2 °° 

<^'')' ~ T(Dn)20n E(2»)'+°-''-'" " ('--)■ (3-2) 

^ I ' n=l 

We evaluate the sum over n in (3.2) by converting it to an integral: 

Af-D/2 foo 

^ ' r(D/2)2^/2 io 
Replacing the sum by an integral is justified when D > 1 because the sum is dominated 
by large values of n lying near the parabola ii? = {D — l)t. Next, we evaluate the integral: 

{R%^^^{2t)'^/' (t^oo) . (3.4) 
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Note that the power of t, which is a critical exponent, for large t is independent of D. 

Although the formula in (3.4) is derived assuming that D > 1, we believe that (3.4) 

remains valid for all -D > 0. We have verified this assertion by means of detailed numerical 

calculations. 

Hausdorff Dimension of a Random Walk 

The Hausdorff dimension D^ of a random walk is defined as the reciprocal of the 



critical exponent of \J {E?)t: the root-mean-square distance from the origin.^ Thus, to 
determine the Hausdorff dimension, we substitute the asymptotic approximation in (3.4) 
with k = 2: 



V{R^)t ~ VDi (t^oo) . (3.5) 

The critical exponent is 1/2 so we conclude that for all dimensions D > we have the 
universal result 

Dn = 2 . (3.6) 

This result is the same as for random walks on Euclidean hypercubic lattices. 
Special Case D = 1 

It is a bit surprising to assert that the dimension of a random walk is 2 even if the 
dimension of the space is less than 2. To demonstrate that this is indeed true, in this 
subsection we evaluate Dh directly for the special case D = 1. When D = 1, we have the 
exact result that for odd t 



4r(| + i) 
v^r(^) 



{R% = l+t+ ^l^^,[ . (3.7) 



Thus, (-R^)t grows like t for large values of t and we have Dh = 2 at D = 1. 
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To derive the result in (3.7) we use an integral representation of the inverse Beta 
function^ to give an integral representation for Cn^t in (2-2): 

4 /"^/2 



Cn,t = - 

n 



r-K/Z 

/ (ix (cosx)* COS [(n — l)x] 
Jo 



with the appropriate factor 1/2 for n = 1. It is important to note that this integral does 
not exhibit {n — t)-parity nor does it vanish when n > t + 1. To compute {R^)t one must 
perform a finite summation of a geometric series. The resulting trigonometric functions 
can then be done explicitly using formulas like^^ 

2 ^/^.._.^tSm[(t + l)x] 



2 f^' 
1 = — / dx [cosxY 

T^ Jo 



sinx 



Note that when D < there is no t-dependence in the root-mean-square as t ^ oo. 
The result for the root-mean-square is divergent for D between and -2; thus Dh = in 
this region. When D < —2, Dh = oo. 
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IV. TEMPORAL MOMENTS OF RANDOM WALKS 

The quantity Hd, as defined in (1.16) represents the probabihty that a random walker 
who begins in the central region will eventually return to the central region. In Ref. 1 it 
was shown that when R^ = n, II d can be computed in closed form as a function of D 
[see (1.18)]. The probability Hd may be regarded as the zeroth temporal moment of the 
discrete probability distribution P2t- 

Another interesting quantity that characterizes the structure of random walks is T, 
the mean time elapsed before the Grst return to the central region. The mean elapsed time 
T is a weighted average of the quantities P2t'- 

oo 

T^^ • (4.1) 

The quantity T is the Grst temporal moment of the probability function P2t divided by 

Although the sum in (1.16) for Uj^ converges for all values of D, it is clear from (2.47) 
and (2.48) that T, as defined in (4.1), converges only if D < or D > 4; T is infinite when 
D lies in the range < -D < 4. Note that when D <2, a, random walker who begins in the 
central region will eventually return to the central region with a probability of unity but 
that the average time for this return is infinite! The reason for this apparently paradoxical 
behavior is that T is an average over all random walks that return to the central region; 
many of these random walks are infinitely long. These long walks dominate the average 
and cause it to diverge. As D increases the shorter walks begin to dominate the first 
moment of P2t', in higher dimension the longer walks almost never find their way back to 
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the central region. Thus, it is not surprising that eventuaUy T becomes a finite quantity. 
The result that T is infinite when < D < 4 is true of Euclidean lattices (see Ref. 4) . 

When L> < it is possible to calculate T in closed form using the reciprocity theorem. 
We multiply (2.37) by 2t and sum over t: 



J2 2^^2t 



t=i 



t=i 


Cl,2t) 


2-D 


= / ,[(2^ ~ 2)Ci,2t-2 — 2tCi,2t + 2Ci,2t-2] 

t=l 


= lim (-2tCi,2t) 


+2ai-D) , 

2-D 



(4.2) 



2-D 



where we have used (1.14). Since 2 — D > 2, then by the asymptotic behavior in (2.26) 
the limit vanishes and we have the exact result 



T = 2C{1-D) {D<0) 



(4.3) 



[Note that we have used the result^ that He) = 1 for D < 2 to obtain (4.3).] 

To calculate T when il> > 4 we use (1.15) to relate the first moment of P2t to the first 
moment of Ci^2t'- 



T 



d 



2t 



t=0 



x=l 



dx Hd 
d -1 

oo 
t=0 



x=l 



,2t 



vt=0 



Ec'i,2tj n 

oo 



t=0 



[({D - 1) - 1]C{D - 1) 
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(4.4) 



where we have used (1.14) and (1.18). 

To evaluate the sum in the numerator of (4.4) we differentiate the continued-fraction 
representation in (1.9). In general, of course, the derivative of a continued- fraction repre- 
sentation is rather complicated. However, for the specific continued fraction in (1.9) we 
find that if we first truncate the continued fraction after the term QnX and then differen- 
tiate once with respect to x, we obtain a sequence of converging structures that we easily 
recognize: 

oo oo ^ oo 

E tC\a< = E ^5^ + E ""-'((O - 1. nf , (4.5) 



71 ■" 
t=0 n=2 n=l 



where 

oo ^ 

k=n+l 

is a Hurwitz function (incomplete zeta function). It is easy to show that the expression 
converges for D > 4 and diverges for D < 4 (it is logarithmically divergent at D = 4) . 
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V. RANDOM WALKS FROM region n TO region m 

Until now we have always considered random walks that begin in the central region 
{region 1). However, to prepare for our consideration of intersecting random walks in 
the next section we address in this section the problem of random walks that originate in 
region n. We will calculate here the probability that a random walker from region n will 
eventually reach region m, where m < n; let H"^'"^ represent this probability. 

We begin our discussion by considering the special case m = 1 and calculate H^'"^, 
the probability that a random walk that begins in region n eventually reaches the central 
region. Note that there are two values of n for which we already know the answer: 

n^'^ = 1 (5.1a) 

and 

U^^^ = Ud , (5.16) 

where Hd is given in (1.18). Equation (5.1a) holds because a random walk that originates in 
region 1 has already reached region 1. To understand (5.1b) recall that Hd represents the 
probability that a random walk that begins in region 1 will eventually return to region 1. 
On the first step of such a walk the random walker goes from region 1 to region 2 with a 
probability of 1 because Pout(l) = 1; now, the probability that the walker reaches region 1 

is ni'2. 

Assume that the random walk begins in region n. After one step, the walker moves 
outward with a probability of Pout(^) to region n + 1 or inward with a probability of Pin{n) 
to region n — 1. Hence the probability that the walker eventually reaches region 1 is a 
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sum of two terms: 

n^'" = Pi„(n)ni'"-i ^ p^^^(^)ni.-+i . (5.2) 

This equation is an ordinary linear second-order homogeneous difference equation whose 
solution is uniquely determined by the two initial conditions in (5.1). The form of the 
exact solution is particularly elementary: 

oo D-D-l 

EfM 

n^'- = t" 'o-^ ■ (5.3) 

Now we consider the general case of a random walk from region n to region m. Again, 
we identify two initial conditions. Clearly, 

n"^'"^ = 1 (5.4a) 

because a random walk starting in region m has already reached region m. Next, we 
calculate n'^'"^+^ following the same procedure used to calculate Hd and we obtain the 
probability in the form of a continued fraction like that in (1.17): 

We then convert this continued fraction to a series like that in (1.17): 

oo tdD-I ^ ' 

Next we construct a difference equation satisfied by n"^'"^. We argue that in the 
first step of a random walk that starts in region n the walker either moves outward with a 
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probability of -Pout(^) to region n + 1 or inward with a probability of Pin{n) to region n — 1. 
Hence the probability of reaching region m eventually is the sum of two probabilities: 

n"^'" = Pin(n)n"^'"-i + Pout(n)n"^'"+i . (5.5) 

The exact unique solution to (5.5) satisfying the two initial conditions in (5.4) is a trivial 
generalization of that in (5.3): 

oo tdD-1 

EfM 



T)D-1 



assuming that m < n. 

Note that if we take Rn = n then when D < 2 the upper end of the sums in (5.6) 
dominate and we have n"^'"^ = 1. When D > 2, n"^'"^ is the ratio of two incomplete zeta 
functions: 

n-» = ,^'^-;-"-^' . (6.7) 

C{D-l,m-l) 
where C(a,n) is defined in (4.6). 
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VI. INTERSECTING RANDOM WALKS 

In rigorous mathematical studies of self-interacting scalar quantum field theory one 
defines a regulated version of the quantum field theory on a Euclidean space-time lattice. 
One then obtains the continuum quantum field theory in the subtle and difficult limit in 
which the lattice spacing tends to zero. In the lattice version of the quantum field theory 
one can show that there is a deep connection between random walks and self-interacting 
scalar quantum field theory (see Ref. 4). Indeed, on the basis of studies of random walks, 
it has been shown that a (j)^^ quantum field theory becomes free when the dimension 
of space-time exceeds ^^i - Proofs that a quantum field theory is noninteracting rely on 
studies of intersecting random walks. ^^ Specifically, it has been shown that a cjP'^ quantum 
field theory is free whenever there is vanishing probability that K independent random 
walks will intersect. 

Rigorous treatments of quantum field theory on a Euclidean lattice have the shortcom- 
ing that the lattice can only be defined in integer dimension. (Thus, while it is known that 
a (j)'^ theory is free in five or more dimensions and interacting in three or less dimensions, 
it is not known if the theory is free in four dimensions.) However, we have shown that on 
a lattice consisting of nested hyperspheres one can define a random walk in any noninteger 
dimension. Thus, if we construct a quantum field theory on such a lattice (as unrealistic 
as such a field theory may be) we may be able to study the transition of this field theory 
from interacting to free as a function of the space-time dimension. The discussion of such a 
rotationally-symmetric quantum field theory is the subject of a future paper. However, in 
this paper we address the crucial question of the probability of intersection of independent 
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random walks. 

The random walks described in this paper are one dimensional; random walks occur in 
the radial coordinate only. Thus, the intersection of two rotationally-symmetric random- 
walk paths has a different meaning from the intersection of two paths on a hypercubic 
lattice. It is intuitively clear, however, that we can study the analog of the intersection of 
two paths provided that appropriate weighting factors are introduced in order to take into 
account the difference between the two kinds of random walks, as we will now argue. 

On a hypercubic lattice one can define^ the conditional probability -P(xi, t, xq, to) for 
a random walker to be on site x at time t starting from the initial position xq at time to- 
Similarly, on a rotationally-symmetric lattice we can define C"'^ as the probability of a 
walker starting at t = in region p to be in region n at time t. Note that C"' is the 
same as Cn,t as defined in Sec. 1. 

Roughly speaking, one can think of the probability Cn,t as a sum over the conditional 

probability P(xi, t,xo, to): 

J2 i^(xi,t,xo,0) ^ 

-XI ^region n ^ ' 

XQ^region 1 

Similarly, one can think of C^'^ as the sum 

J2 i'(xi,t,xo,0) ^ 

x-j^ E^region n ^ '^ 

XQ^region p 

The probabilities C"'^ are not mutually exclusive. Thus, the sum J2t^p-n ^t'^ ^ where 
we take n < p, is not itself a probability; rather, it represents the mean time spent in the 
nth region. This sum is finite only for D > 2 [see, for example, Eq. (1.14)]. One can 
define^ -Pi(xi, t, xq, to) as the probability for a walker on a hypercubic lattice to be at site 
X at time t for the first time starting from the initial position xq at time to. Similarly, 
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P"'^ represents the probability of a random walker starting in region p to be in region n 
{n < p) at time t for the first time. Roughly speaking, P[^'^ can be thought of as the sum 

J2 A(xi,t,xo,0) ^ 

x.-\^ ^region n ^'^ 

XQ^region p 

The mutually exclusive probabilities P2t defined in (1.15) are given by P2t = -P2t-i- 
In contrast to J2t^p-n ^f'^? the sum on Pj^'^ is a probability and in fact H'^'P (n < p), as 
defined in Sec. V, is given by 

oo 

J2 Pt'" = n"'P . (6.4) 

t^p—n 

The generalization of the generating-function formula (1.15) is found to be: 

cx) ^oo ^n,p t 

2_^ ^t ^ - ^oo ^n,n . • 10.5J 

To obtain (1.15) from (6.5) note that C\'^^ = C^^ 

On a hypercubic lattice, the probability P(xi,X2,A) that two random walk paths 
starting at xi and X2 (xi ^ X2) will intersect in a region A, satisfies^ 



p(xi, X2, A) < 5Z n(z - xi)n(z - X2) 



(6.6) 

z6A 

where n(x) = Xl^o ^i(^' ^)- 

One is interested in the continuum limit of this expression. In order to calculate 
the right side of (6.6) we use our results for the rotationally-symmetric random walk 
probabilities H"^'" in (5.6) and we note the difference between n(z — x) and 11'^'^, which 
includes the summation over the sites in the regions. The appropriate probability that 
must be used in (6.6) is 11'^'^ multiplied by 
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The right side (6.6) is thus given by: 

j^n-n-W^i^ • (6-7) 

To obtain the continuum hmit we now rescale aU dimensions by 

p -^ ap, q -^ aq, N -^ aN . (6-8) 

From (5.6) one notes that H^'P ~ (p/n)'^~^ and the correction factor in (6.7) ehminates 
the n dependence for large n in the probabihties entering into the right side of (6.6). After 
rescahng (6.8) we find that (6.7) is given by 

/oiN 
rfnn^-in(2D-^)n2(2-^) , (6.9) 

where we show the exphcit n dependence of the different factors in (6.7). Thus, the right 
side of the inequahty (6.6) vanishes for large a as a^~^. 

More generally, for the case of K independent random walks we find that the right 
side of (6.6) is proportional to 

/aN 
rfn n^-1 = a(2-^)^+^ • (6.10) 

Thus, the probability that K random walk paths intersect vanishes if 

2K 
D>^^ . (6.11) 

K-1 ^ ' 

This reproduces the well-known condition for a self-interacting scalar quantum (/)^^ field 
theory to be nontrivial. It gives the critical dimension D — A for (f)^ and D — ?> for (f)^ field 
theories. 
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FIGURE CAPTIONS 

Figure 1. Three-dimensional plot of Cn,t on an (n, t) lattice for the case D = 3. For a given 
value of n > 1, n — 1 is the smallest value of t for which Cn,t is nonzero. (This is because 
the random walker cannot reach region n until time n — 1.) We find that when D > 1, 
Cn,t at first increases with t, reaches a maximum, and then dies off to as t continues 
to increase. Also, for fixed t as we vary n, we again observe a local maximum. Thus, in 
the plot of probabilities there is a crest; this crest lies on a parabola: t{D — 1) = n^. The 
parabola is indicated on the graph. The graph shows that the crest is very broad. 

Figure 2. Same as in Fig. 1 except that D = 1/2. When < D < 1 the maximum value 
of Cn^t for each fixed value of t lies along the t axis; it no longer lies on a parabolic curve. 
As t ^ oo with n fixed, Cn,t tends to zero very slowly. 

Figure 3. A plot of Ci t (solid curve) and the leading asymptotic approximation to Ci^t in 
(2.27) (dashed curve) for t = 2, 4, 6, 8, . . . , 3000 for the case D = 1/2. The asymptotic 
approximation clearly becomes more accurate as t increases. 

Figure 4. Same as in Fig. 1 except that D — —1. When D < we observe a new behavior 
not present in Figs. 1 and 2; namely, that as t ^ oo with n fixed, Cn^t approaches Cn, oo, 
a nonzero time-independent function of n. The function Cn,oo is given in (2.36). 
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